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ANALYTICAL  SOLUTION  OF  ELASTIC-PLASTIC  THICK-WALLED 
CYLINDERS  WITH  GENERAL  HARDENING 


Peter  C.T.  Chen 

U.S.  Army  Armament  Research,  Development,  and  Engineering  Center 
Close  Combat  Armaments  Center 
Benet  Laboratories 
Water vl iet ,  NY  12189-4050 


ABSTRACT .  This  paper  presents  an  analytical  solution  for  the  elastic- 
plastic  behavior  of  thick-walled  cylinders  under  internal  pressure.  The  general 
hardening  law  employed  in  this  investigation  is  a  piecewise  linear  represen¬ 
tation  of  arbitrary  stress-strain  curves  in  uniaxial  form.  Closed-form  analyti¬ 
cal  solutions  are  developed  for  the  stresses,  the  elastic  and  plastic  strains, 
and  the  displacements  by  using  Treses' s  yield  criterion  and  its  associated  flow 
rule.  ExDerimental  data  obtained  from  cylinders  made  of  either  SAE  1045  steel, 
OFHC  copper,  or  aluminum  alloy  1100  are  used  to  determine  the  material 
constants.  Numerical  results  for  partially-plastic  and  fully-plastic  cylinders 
are  presented  for  the  radial  distributions  of  plastic  hoop  strain,  radia1 ,  and 
tangential  stresses. 

1.  INTRODUCTION.  Of  all  the  available  elastic-plastic  solutions,  the 
Droblem  of  pressurized  thick-walled  cylinders  has  received  the  greatest  atten¬ 
tion.  This  is  because  of  the  symmetric  nature  of  the  problem  and  its  practical 
imDortance  to  pressure  vessels  and  the  autofrettage  process  of  gun  tubes.  Many 
solutions  for  this  problem  have  been  reported  over  the  last  four  decades  [1-3]. 
Analytical  solutions  can  be  obtained  only  when  simplifying  assumptions  are  made 
regarding  material  properties.  Bland  [2]  developed  analytical  solutions  for 
materials  with  linear  hardening  proDerties.  Recently,  Megahed  [3]  considered  a 
nonlinear  hardening  law  a  -  Y  +  A*eDn  in  uniaxial  form  and  developed  an  approxi¬ 
mate  solution  for  any  value  of  the  strain-hardening  exponent  n.  Closed-form 
analytical  solutions  for  the  plastic  hoop  strain  can  be  obtained  only  for  four 
particular  values  (n  *  1,  1/2,  1/3,  and  1/4),  and  the  integral  has  to  be  eval¬ 
uated  numerically  for  n  =  1/3  and  1/4. 

The  general  hardening  law  employed  in  this  investigation  is  a  Diecewise 
linear  representation  of  actual  stress-strain  curves  in  uniaxial  form.  A  finite 
number  of  straight  lines  can  represent  arbitrary  curves  with  greater  accuracy 
than  other  representations  [4],  The  problem  is  formulated  in  a  manner  similar 
to  [2,3]  by  using  Tresca's  yield  criterion  and  the  associated  flow  rule. 
Closed-form  analytical  solutions  are  developed  for  the  stresses,  the  elastic  and 
plastic  strains,  and  the  displacements. 

2.  BASIC  EQUATIONS.  Consider  a  long  thick-walled  cylinder,  internal 
radius  a  and  external  radius  b,  that  is  subjected  to  internal  pressure  o 
causing  partial  plastif ication.  Assuming  small  strain  and  no  body  forces  in  the 
axisymmetric  state  of  generalized  plane-strain,  the  radial  and  tangential 
stress,  <7r  and  <jg,  must  satisfy  the  equilibrium  eauation 

r(daP/dr)  -  aQ  -  ar  (1) 

and  the  corresponding  strains,  er  and  eg,  are  given  in  terms  of  the  radial 
displacement,  u,  by 


€r  *  du/dr  ,  eg  =  u/r 


(2) 
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Total  strains  are  decomposed  into  elastic  and  plastic  components  and  the  strain- 
stress  relations  are 


er  =  far  ~ 

v(a0+<iz)]/E  +  €rP 

(3a) 

1 

<X> 

£ 

tl 

<t> 

U) 

V(Or+Oz) ]/E  +  £gP 

(3b) 

0\ 

N 

II 

i — i 

Q 

N 

1 

V(dr+<7g)]/E  +  £ZP 

(3c) 

where  E  and  v  are  elastic  constants. 

Subject  to  a q  Z  oz  }  or, 

Tresca's  cri- 

terion  states  that  yielding  occurs  when 

<tq  -  ar  =  a(£P)  (4) 

where  the  yield  stress  a  is  a  function  of  plastic  strain  ep.  The  associated 
flow  rule  states  that 


Hence,  from  Eq.  (3c) 


d£gP  =  -d£rp  *  0 


and  dezp  =  0 


(5) 


O 2  =  u(cr+cg)  +  Eez 

and  the  total  axial  force  on  any  section  is 

F  =  2nva2p  +  7rE(b2-a2)£z 


(6a) 


(6b) 


There  are  three  cases  of  importance:  first,  plane-strain,  ez  =  0;  second,  a  tube 
with  open  ends,  F  =  0;  and  third,  a  tube  with  closed  ends,  F  =  jra2p.  In  the 
latter  two  cases,  substitution  into  Eq.  (6b)  determines  ez.  Since  ez  is  now 
known,  Eos.  (3a)  and  (3b)  are  inverted  in  order  to  express  stresses  in  terms  of 
total  strains  and  plastic  hoop  strain  as  follows: 


12)  yields  the  following  differential  equation: 


d2u  1  du  u 
dr2  +  r  dr  "  r2 


Integrating  with  respect  to  r  leads  to 


1-2V  2£g' 


du/dr  ♦ 


where 


J  =  /  £g  r-1dr 
a 


♦  ”e2] 

(7a) 

+  vcz] 

(7b) 

and  (7b)  into  Eqs. 

( 1 )  and 

~dr“J 

(8) 

+  2C 

(9a) 

(9b) 
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Integrating  again  yields  the  analytical  solution 


u  =  -(i---)rJ  +  Cr  +  D/r 
1-v 


(10) 


where  C  and  D  are  integration  constants  to  be  determined  from  boundary  con¬ 
ditions;  Up  *  -d  at  r  =  a  and  <rr  =  0  at  r  =  b. 

Uoon  substitution  of  the  resulting  values  of  C  and  D  into  the  expressions 
of  displacement ,  radial  and  hoop  stresses,  the  following  distributions  are 
obtained: 

ar  =  -p  -  EJ  +  (P+EJ0) (l-a2/r2)/(l-a2/b2)  (11a) 

c70  =  -p  -  E(J+e0P)  ♦  (p+EJ0) (l+a2/r2)/ (l-a2/b2)  (lib) 

EU/r  =  (1+v) (l-2v)or  +  2(l-v2 ) (a2/r2 ) (p+EJ0)/(l  -  57)  -  i>Eez  (11c) 


where  E  =  E/(l-v2),  and  JQ  is  the  value  of  the  integral  J  at  the  plastic  front, 
r  =  0,  i .e. , 

P  P/a  _ 

J0  =  /  e0p  r-i  dr  =  /  eeP  €_1d«  (lid) 

Note  that  £0°  =  0  and  J  =  JQ  throughout  the  outer  elastic  zone_defined  by 
p  <  r  <  b.  At  the  plastic  front,  the  Tresca  effective  stress  a  -  ctq  -  ar  =  Y, 
where  Y  is  the  initial  yield  stress  and  also  £0P  *  0.  Therefore,  using  Ea. 

(lib)  to  provide  a,  one  readily  obtains 


-  Oq  pi  a2 

p  ♦  E  J0  =  (1  - 


(12) 


Using  Eq.  (7),  the  distributions  of  ar,  oq,  and  u  can  be  written  in  simpler 
forms  as  follows: 


<rr  =  -0  -  EJ  +  |  (p2/a2  -  p2/r2 )  (13a) 

(Jo  =  or  +  Y(p2/r2  -  --  £0P)  (13b) 

°o 

Eu/r  =  (1+v) (l-2v)or  +  Y ( 1-v2 )p2/r2  -  vEez  (13c) 

It  is  obvious  from  Eq.  (13b)  that  Tresca  effective  stress  a  is  simply  given  oy 

a  =  Y (p2/r2  -  f-  £0P)  (1A) 

Therefore,  if  the  radial  variation  of  plastic  hoop  strain  is  known,  the  integral 
J  and  all  field  quantities  can  be  determined. 
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3.  GENERAL  HARDENING.  The  general  hardening  law  employed  in  this  investi¬ 
gation  is  a  piecewise  linear  representation.  Arbitrary  stress-strain  curves  in 
uniaxial  form  can  be  approximated  by  a  finite  number  of  straight  lines  [4].  The 
straight  line  through  the  origin  is  given  by  the  relation 

a  =  Ei  (15) 

where  E  is  Young's  modulus.  All  of  the  other  straight  lines  are  given  by  the 
relation 


a  -  (1-m.j  )o0t  +  m-jEe 


(16) 


where  a0j  is  the  stress  at  the  intersection  of  the  two  straight  lines  given  by 
Eqs.  (15)  and  (16),  and  m-jE  is  the  slope  of  the  straight  lines  given  by  Eq.  (16). 
Let  Oj ,  £j  be  the  stress  and  strain  at  the  intersection  of  two  straight  lines 
with  slope  m-j_iE  and  m-jE  as  shown  in  Figure  la.  Then 

cr-j  =  ( 1— m-j  ^  +  m-j_}Ee-j  —  (1— m^)OQ-j  +  m-jEc-j 

which  leads  to  e-j  and  cr-j  in  terms  of  cr0-j  and  m-j 

Ee-j  =  f  ( 1-m-j  )a0i  -  (l-m.j>1)a0.j_1]/(m.j_i-m.j )  (17a) 

and 

a-j  =  [m-j_i(l-m-j  )a0-j  -  m-j  (l-m-j-!  )<70i-l]/(mi-l-mi )  (17b) 

Eq.  (16)  can  be  written  also  as  a  function  of  effective  plastic  strain  ep  as 
shown  in  Figure  lb. 


a  =  o0-j  +  h-jEe°  e-\°  i  eP  ^  ei  +  l°  (18) 

where  h-j  =  m-j/d-m-j). 

Since  eg13  =  -€r°  and  e2p  =  0,  the  effective  plastic  strain  ep  is  determined 
as  2£qP//3,  and  hence,  Eq.  (14)  is  rewritten  in  terms  of  the  plastic  strain  in 
the  tube  as 


o  =  Y(p*/r* 


-!  I  ip> 

2  Y  1 


(19) 


A  comparison  between  the  expressions  for  effective  stresses  provided  by  Eos. 
(18)  and  (19)  yields  the  following  explicit  equation  for  eP: 


§  JP 

Y 


(p*/r«- C-jJ/b-j 

which  is  valid  in  e-j“  4  eH  4  e ^ ^ ^  and  z  r  }  r-j  +  j  and 

bi  ■  (✓3/2)/(l-v*-)  +  hi  ,  C-j  =  a01/£T0 


(20) 


(21) 
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The  values  of  r-j  and  r^+1  can  be  determined  by 

r i  =  P C b ■?  y  +  C-j]  ^ 


ri+l  s  P[bi  9  ei+l^  +  ^il-J* 


If  £}P  =  0  and  <70i  *  Y,  then  =  p.  This  is  true  for  most  materials.  Since 
r-j  J  a,  the  calculation  of  r-j  for  i  =  1,2... m  should  stop  when  the  above  rela¬ 
tion  is  violated,  i.e.,  rm+1  <  a.  Let  us  define  V-j  for  i  =  1,2,... m  by  the 
following  i ntegr a 1 : 

v,-  *  fr'  Mpr  =  5-C  &-^r- 

ri+l  Y  r  ri+l  r  r 

=  55”  H =?-->*  '  (;~)2  ”  2Ci  ln(ri/ri+i)]  (2 

*Di  ri+l  ri 


f  E  -  dr  t  m  E  -p  dr  ..  ..  rr  E  -o  dr 

V  ”  K  Y  £P  r”  =  /  Y  £  r“  Vm-1  +  - *  ’  Vi+1  +  /  3  C'  r~ 

a  a  r.j+1 


2k  C(i)Z  •  {k]'  - 2C™ *"  r]  +  V-1  ♦•••♦  V-1 


+  25”  -  2C,  in  (r/ri+1: 


and  the  maximum  value  of  V  is 


*0  ■  u  5p  r  -  ik  t«I»!  -  -  2C»  >"  r1  *  i  ^ 


The  integrals  ( i*l ,2, . . . ,m-l) ,  V,  and  V0  given  analytically  by  Eqs.  (23), 
(24),  and  (25)  can  be  easily  evaluated.  The  integral  J  is  related  to  the 
integral  V  by 


J  =  r  e  v 


All  field  quantities  u,  er,  eg,  ar,  eg,  az,  and  eg  can  now  be  calculated. 
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4.  MATERIAL  PROPERTIES.  Test  members  were  made  froe  three  different 
metals  as  follows:  SAE  1045  steel,  OFHC  copper,  and  aluminum  alloy  1100  [4]. 

The  values  of  the  elastic  constants  {E  and  v)  for  the  three  metals  are  shown  in 
Table  1.  The  values  of  the  constants  (aQ.j.  m,  ,  e^)  approximating  the  plastic 
portion  of  the  stress-strain  diagram  for  three  metals  are  shown  in  Table  2. 

TABLE  1.  ELASTIC  CONSTANTS  FOR  THREE  METALS 


Material 

E,  Ksi 

v 

SAE  1045  steel 

30,000 

0.29 

OFHC  copper 

16,000 

0.35 

aluminum  allov  1100 

10,250 

0.33 

TABLE  2.  PLASTIC  CONSTANTS  FOR  THREE  METALS 


Straight  Line 

— 

cr0i  <  Ksi 

— 

mi 

,  Ksi 

e-j ,  % 

SAE 

1045  Steel 

1 

43.4 

0.05083 

43.4 

0.145 

2 

54.0 

0.02858 

66.924 

1.687 

3 

80.0 

0.00847 

90.638 

4.453 

4 

95.0 

0.00309 

103.542 

9.532 

5 

111.0 

0.00128 

122.280 

29.745 

OFHC  Cooper 

1 

2.50 

0.17125 

2.50 

0.016 

2 

3.25 

0.07063 

3.686 

0.059 

3 

4.00 

0.03125 

4.553 

0.136 

4 

5.37 

0.01991 

o 

o 

t 

0.765 

5 

8.40 

0.01313 

14.151 

2.790 

6 

21.0 

0.00450 

27.484 

9.137 

7 

39.0 

0.00078 

42.757 

30.350 

Aluminum  Allov  1100 

1 

8.0 

0.67024 

8.0 

0.078 

2 

11.0 

0.32683 

11.942 

0.135 

3 

13.0 

0.09561 

13.557 

0.184 

4 

14.7 

0.01590 

15.007 

0.332 

5 

16.1 

0.00210 

16.310 

1.131 

Each  of  the  stress-strain  diagrams  can  be  approximated  by  a  finite  number  of 
straight  lines  with  extreme  accuracy.  The  error  introduced  by  the  approximation 
is  less  than  1  percent  for  all  cases. 
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5.  NUMERICAL  RESULTS.  Typical  results  for  the  analytical  solution  are 
presented  first  by  means  of  prescribing  a  plastic  front  and  determining  the 
corresponding  plastic  hoop  strain  and  radial  and  hoop  stresses  in  the  tube.  A 
tube  with  b/a  *  2  is  employed,  and  the  plastic  front  is  prescribed  at  p/a  =  1.0, 
1.2,  1.4,  1.6,  1.8,  and  2.0.  Figure  2  illustrates  the  stresses  and  plastic  hoop 
stra.ns  obtained  using  the  material  constants  for  SAE  1045  steel.  Figures  3  and 
4  present  similar  results  for  OFHC  copper  and  aluminum  alloy  1100,  respectively. 
Figure  5  shows  a  comparison  of  stresses  and  clastic  hoop  strains  for  three 
partially-plastic  tubes  at  p/a  *  1.6.  Figure  6  presents  a  similar  comparison 
for  three  fully-plastic  tubes  at  p/a  »  2.0.  Future  work  related  to  the  results 
obtained  here  will  look  into  the  elastic-plastic  behavior  of  the  tube  during 
pressure  release.  The  influence  of  phenomena  such  as  the  Bauschinger  effect  on 
residual  stresses  should  be  modelled  [5,6]. 
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